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Abstract 

Using the Ehlers transformation along with the gravitoelectromagnetic approach to stationary 
spacetimes we start from the Morgan-Morgan disk spacetime (without radial pressure) as the seed 
metric and find its corresponding stationary spacetime. As expected from the Ehlers transformation 
the stationary spacetime obtained suffers from a NUT-type singularity and the new parameter 
introduced in the stationary case could be interpreted as the gravitomagnetic monopole charge (or 
the NUT factor). As a consequence of this singularity there are closed timelike curves (CTCs) in 
the singular region of the spacetime. Some of the properties of this spacetime including its particle 
velocity distribution, gravitational redshift, stability and energy conditions are discussed. 
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I. INTRODUCTION 



Finding a new solution of Einstein field equations from a given solution of the same field 
equations has a long history and mathematically is related to the problem of finding a solu- 
tion to a set of nonlinear partial differential equations from another known solution. Here we 
make use of a very special transformation introduced by Ehlers [1] through which one could 
start from a static spacetime and end up with a stationary solution. We apply this trans- 
formation to two static spacetimes, namely Schwarzschild metric and Morgan-Morgan disk 
space [2] and in the meantime employ the concepts introduced in the gravitoelectromagnetic 
approach (1+3 formalism) to spacetime decomposition to yield a physical interpretation of 
the spacetime found through the transformation. The use of the gravitoelectromagnetic con- 
cepts are justified by the fact that the main equation of the Ehlers transformation includes 
the gravitomagnetic potential of the stationary spacetime. The outline of the paper is as 
follows. In section II we introduce briefly gravitoelectromagnetism and the Ehlers transfor- 
mation and use the latter to derive, as an example, the NUT solution [3] from Schwarzschild 
spacetime and also Chazy-Curzon-NUT (CC-NUT) solution from Chazy-Curzom static met- 
ric [4,5]. In section III we introduce briefly the static Morgan- Morgan disk space (without 
radial pressure) and some of its properties. In section IV using Ehlers transformation we 
find the corresponding stationary spacetime representing the spacetime of a finite disk of 
mass M and gravitomagnetic charge I. In section V we discuss some of the properties of this 
stationary disk including its particle velocity distribution, gravitational redshift, stability 
and energy conditions. In the last section we summarize the results obtained. 

II. GRAVITOELECTROMAGNETISM AND EHLERS TRANSFORMATION 

The (1 + 3)-decomposition (threading) of a spacetime by a congruence of timelike curves 
(observer worldlines) leads to the following splitting of the spacetime interval element [6]; 



where dL and dT are defined to be the invariant spatial and temporal length elements of 
two nearby events respectively. They are constructed from the normalized tangent vector 
u a = 4t to the timelike curves in the following way j^j; 



ds 2 = dT 2 - dL 2 , 



(1) 



dL 2 = h ab dx a dx b , 



(2) 
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dT = u a dx a , (3) 

where 

hab = -g a b + u a u bl (4) 

is called the projection tensor. Taking h = |£| 2 and A a == — equations (1) and (4) can 
be written in the following alternative forms; 

ds 2 = h(A a dx a ) 2 - h ab dx a dx b ; h ab = -g ab + hA a A b . (5) 

Using the preferred coordinate system in which the timelike curves are parameterized by the 
coordinate time x° of the comoving observers; 

C = (1,0, 0,0) ; A a = (-1,-— ), (6) 
and the above spatial and spacetime distance elements will take the following forms [6,7]; 

dL 2 = dl 2 = laP dx a dx^ (7) 

ds 2 = e 2U (dx° -A a dx a ) 2 -dl 2 , (8) 

where 

e 2U = goo , A a = -gojgoo, (9) 

and 

/ . goagop-, / in x 
goo 

Introducing U and A as the gravitoelectric and gravitomagnetic potentials respectively, their 
corresponding fields are j3] [7]; 

E s = -W, (11) 

B g = curlA, (12) 

where the differential operators are defined in the 7 space. On the other hand Ehlers 
transformation, also called gravitational duality rotation, in its original form states that if 

g mn dx m dx n = e 2U {dx ) 2 - e~ 2U dl 2 , (13) 

with dl 2 = e 2U dl 2 representing the conformal spatial distance, is the metric of a static 
exterior spacetime, then, 

g mn dx m dx n = (acosh(2U)y 1 (dx° - A p dx p ) 2 - acosh(2U)dl 2 , (14) 



(with a = constant > 0, U = U(x a ) and A$ = Ap(x a )) would be the metric of a stationary 
exterior spacetime provided that A a satisfies the following equation; 



aVie af3v U> v = A [a ^, (15) 

where 7 = det7 Mi ,, 7^ being the conformal spatial metric. Hereafter we call the above 
equation the Ehlers equation. 

Considering the Ehlers transformation and in particular Ehlers equation (15) in the context 
of gravitomagnetism, then finding a stationary solution corresponding to a given static 
solution, amounts to finding the gravitomagnetic potential A a related to a given static 
(potential) function U. In what follows, as an example of the procedure, we derive the 
NUT space [8] as the stationary spacetime with a radial gravitomagnetic field which is the 
gravitational dual of Schwarzschild space. 

A. NUT space from Schwarzschild through Ehlers transformation 

Starting from the schwarzschild mertic in the form (13) in which ; 

U c = - ln(l - 2m/ r) + C, (16) 
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and 



dl 2 = dr 2 + r 2 exp{2U){d6 2 + sin 2 6d<f) 2 ), (17) 



we find that in this case the Ehlers equation (15) reduces to, 

ar 2 e 2U sm6'y rr U r = ~A <l>> e, (18) 

in which we set Ag^ = to keep the resulted spacetime single valued and axially symmetric. 
Using (16) the gravitomagnetic field is found to be, 

A^O) = 2am cos6. (19) 

So the resulted stationary metric is given by (14) as; 

ds 2 = r(r ~ 2m) (dt - 2am cos Bd^f - dr 2 - af(r)dn 2 , (20) 

afyr) r[r — 2m) 



where 

f(r) = r 2 (- 

2ci 



1 + c 2 

f(r) = r 2 ( -) + 2m 2 d - 2mrc x , c x = e 2C . (21) 



This is the well known NUT space [3,11]. To recover the Schwarzschild spacetime metric as 
a^Owe need 

«=j^ , Ici|<l- (22) 
Now applying the following changes of variables and redefinitions of constants ; 

1 - c 2 

M = m(- \) , am= I , (r - la) = R, (23) 

1 -\~ C-y 

to the line element (20) we have, 

^ = R2 ~™ + R i l \ dt - 2/cos^0) 2 - R2 * 2 2 MR-P dR2 - ^ + /2) ^' (24) 

which is the NUT space in its more common form [8,11,12]. Note that having started from 
a potential with C = (or a = 1) we would have ended up with the pure (M = 0) NUT 
space. In other words although adding the constant C in the static case does not correspond 
to a different spacetime, in the Ehlers transformation it does. This actually answers another 
question of concern, namely, what is the seed metric (in Ehlers transformation) of pure NUT 
space?. By the above argument for both NUT and pure NUT spaces the seed metic is the 
Schwarzschild metric. 



B. stationary spacetimes with A a = A a (r,9) 

As another example of the application of the Ehlers transforemation in this section we 
consider a family of solutions for a gravitomagnetic potential of the form 

A = A M)0, (25) 

(In the spherical or the schwarzschild-type coordinates (t, r, 9, 0) ) which according to the 
1 + 3 formalism [7-9] could, in principle, produce stationary spacetimes with gravitomagnetic 
field components B g (r, 8)r and B g (r, 9)9. Using Ehlers transformation one could easily show 
that for an A of the type (25) one ends up with the equation V 2 U — for U with the 
general solution [11] 

U = Z n=0 a n r-( n+ Vp n (cos9). (26) 

For the simplest case, i.e. n = we have U oc - which, in the static case corresponds to 
the Chazy-Curzon metric [4,5] and in the stationary case (through Ehlers transformation) 
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to a spactime dual to Chazy-Curzon with a gravitomagnetic potential of the form A^iO) = 
acosO + b [13,14]. It is notable that this has exactly the same form as the gravitomagnetic 
potential of NUT space and consequently the parameter a could be interpreted as the NUT 
factor (or magnetic mass). Another example in which Ehlers transformation is examined 
on a static spacetime to yield an stationary space is the case of cylindrical NUT space [15] 
where it is shown that it could be obtained from Levi-Civita's static cylindrically symmetric 
solution. 



III. STATIC SPACETIME OF A PRESSURELESS (MORGAN-MORGAN) DISK 

Static axially symmetric spacetimes in the Weyl canonical coordinates are given by the 
following general form [11]; 

ds 2 = e~ 2U [e 2k {dp 2 + dz 2 ) + p 2 ^ 2 ] - e 2U dt 2 , (33) 

in which U and k are functions of p and z. Morgan and Morgan [2,11] found the following 
solution for a finite disk (of mass M and radius a) without radial pressure in terms of the 
oblate ellipsoidal coordinates (£,77); 



U = (arccott + ^[(3£ 2 + l)arccot£ - 3£](3r/ 2 - 1)^ 



(34) 



and 



k = - A M 2 p 2 a~ 4 \(^) 2 B 2 (0 - (1 + V 2 )A 2 (0 - 2£(1 - r?)A(Z)B{zj\ , (35) 

A(0 = £arccot£ - 1 , 5(0 = ^[^^ - axccotf], (36) 
in which the connection between the ellipsoidal (£, rf) and Weyl (p, z) coordinates; 

p 2 = a 2 (l + £ 2 )(l-r/ 2 ) , z = atrt, \r,\ < 1 , < £ < 00, (37) 

infer that the disk is located at £ = , \rj\ < 1 . In [16] the authors have found the form of 
the same metric in Weyl canonical coordinates as the imaginary part of a complexified bar 
with the functions U and k as follows; 

M , .Re\R] - ia. 

U = , 38 

2ia l Re[R]+ia h v ; 



where 
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1M ftl^ 

2 2 R K ' 

where R = p 2 + (z — ia) 2 . Exact thin disk solutions of Einstein field equations are dis- 
cussed extensively in different contexts in the literature. The above static solution, first 
derived by Morgan and Morgan [2], is considered to represent the spacetime of a finite disk 
of counterrotating particles. Properties of static counterrotating disks have been studied in 
[17,18]. Superposition of a central black hole with a surrounding thin disk of counterrotat- 
ing particles is discussed in [19] and infinite self-similar counterrotating disks were studied 
in [20,21]. In the following section using the Ehlers transformation we obtain the station- 
ary spacetime corresponding to the Morgan-Morgan static disk space. To our knowledge 
the exact NUT extension of the above Morgan-Morgan disk space has not appeared in the 
literature. 



IV. STATIONARY SPACETIME OF A FINITE NUTTY DISK 

Starting with the Morgan-Morgan static disk space given by equations (34-36) and using 
the fact that 3r] 2 — 1 = 2P 2 (r]) we write (34) as follows 

U c = ~ ^arccot£ + i[(3£ 2 + l)arccot£ - 3£]P 2 (??)^ + C. (40) 

Now choosing an axially symmetric gravitomagnetic potential A = A^^r))^, by Ehlers 
transformation we have; 

A^ = 2aa{l + e)U^ (41) 

A u = -2cm(l - rj 2 )U v , (42) 

where a is the Ehlers transformation parameter (duality rotation parameter) and a is the 
radius of the disk with mass M. Since the right hand sides of the above two equations 
include Legendre polynomials Pq{t}) and P 2 {r]) we can expand the gravitomagnetic potential 
as follows; 

Mr,0) = Mv)Pi(0 + fs(v)Ps(0- (43) 
By substituting the above potential into the equations (41) and (42) we end up with the 
following solution for the gravitomagnetic potential; 

MZ> V) = 3^(1 - ^)(1 + e 2 ) arctan£+ 

\rfi [K 3 - 2e - 4/3 + ir£)v 2 - - 2^ + + d, (44) 
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where as before I = Ma is the NUT factor. To regain the Morgan-Morgan space as I — > 0, 
we choose a = 2c\ where c\ = e 2C , so that the line element of the stationary spacetime of a 
finite disk with mass M and magnetic mass (NUT factor) I is given by ; 

d$2 = { ^ e2U + e_2t/)-1 [dt ~ M ^ V)d<l)]2 

~ + e~ 2U ) + V 2 )(^ + YT^f) + (1 + e)(l - . (45) 

This is verified to be an exact solution to the vacuum Einstein field equations using the Maple 
tensor package J3j. As the original Morgan-Morgan disk space is asymptotically flat, it is 
not difficult to see that the above stationary spacetime shares the same property. Hereafter 
we call the above solution the MM-NUT disk space. The presence of the Dirac-type (NUT- 
type) singularity could be inferred from the form of the gravitomagnetic potential (44), as 
one could not make it vanish simultaneously for both the positive (z > = rj = 1) and 
negative (z > = rj = —1) half-axes (p = 0) by any choice of the constant C%. In other 
words only half of the axis could be made regular, either the positive half-axis (for C\ = 21) 
or the negative half-axis (for C\ = —21). Therefore a shrinking loop around the singular 
half- axis will have a non-zero circumference even when p — > 0. Now consider a t — Constant 
hypersurface and take the lower half-axis (n = —1) to be singular ( i.e. C% = 21). In this 
case the coefficient of dcfi 2 in (45) becomes positive at sufficiently small p (n w —1); hence <ft 
becomes a timelike coordinate and being a cyclic variable with a period of 2tc, this means that 
the spacetime contains closed timelike curves in the singular section; another characteristic 
of NUT-type singularity. 

V. PARTICLE VELOCITY, STABILITY, ENERGY CONDITIONS AND GRAVI- 
TATIONAL REDSHIFT 

In this section we calculate particle velocity distribution of the MM-NUT (stationary) 
disk found in the previous section. We also discuss its stability and energy conditions and 
obtain the gravitational redshift suffered by the photons emitted from different radii on the 
disk. Hereafter we take a disk with unit radius i.e. we set a = 1 for convenience. 
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A. Particle velocity 



The general form of the metric for stationary axially symmetric spacetimes (in cylindrical 
coordinates) is given by; 

ds 2 = g tt (p,z)dt 2 + 2g t4> (p, z)dtdcj) + g aa (p, z)(dx a ) 2 a = p,z,<p. (46) 

For a zero angular momentum observer (ZAMO) in these spacetimes, the velocity of a 
particle in a circular motion (in the p = Constant, z = plane), is given by [22]; 



d<f> . 9t4> 



2 



vl= 2 ^ (77 + — ■ (47) 

Therefore we need to find -i, for which, we use the energy-momentum tensor of the particle 

distribution which in turn could be calculated by the formalism of distributions [19,23]. 
Using the general form of the energy-momentum tensor of a system of particles [7], 

TT = Vmtl-gyi^^Six* - x), (48) 

and the fact that ^ = (-^) 1//2 we end up with the folowing final result ; 

n2 _ T+* _ b pp {gPPgW) + b tt (g tt g^ - g^g*) (4g) 



v dt ' T a b pp (g tt gpp) + b H (g u g^ - g^gW) 

in which, through the Einstein field equations, the components of the energy-momentum 
tensor are given by j^ . 

V = ~ \b az 5l - b zz 5 a b + g az 5 z b - g zz b a b + b c c (g zz 5 a b - g az 5 z b )} 6(z), (50) 

where b ab are the components of the discontinuity in the first derivative of the metric tensor 
(at z = 0), 

b a b — 9ab,z\z=0+ ~ 9ab,z\z=0~ = ^9ab,z\z=0+- (51) 

Computing these elements for the metric (45) we have; 



b tt = 30MVl-^-(p) 



(53) 



b H = 3QMy/T=?F-(p) (p 2 + |^|) , (54) 



where 

F±(p) = ^f 2 e^ (p2 - 2) ± e-f ^- 2 ), (55) 

and 

G(p)=4(^ = 0) = (2/-2/[l-p 2 ]f) 2 . (56) 

Now substituting (52-56) into (50) and use the outcome in (49) and (47) we obtain the 
follwoing expression for the particle velocity, 

*♦ = WM (Fl[fi)pl ~ G(P)) + W-Uo ) • (57) 

where 

. . i— — ( I 4 2 / 3Mvr(p 2 -2) x t I 4 / 3Mtt(p 2 -2) x 
fl(p) = 4^^— p 2 exp( ^ ') +4 — exp( ^ ^ 

72 

+ 24vr— p 2 - 192lnG?(p) + 5767r/ 2 p 2 (l - p 2 ) + 1927r/ 2 p 6 

„ n 2 , 3Mtt(p 2 -2V 64 . 3Mvr(p 2 -2)A~* 
+ 48vrp 2 exp( ^ >-) - — exp( ^ ^ J . (58) 

For the disk to be physical we impose the condition that the velocity of the particles at the 
rim of the disk do not exceed that of the light i.e., 

v*(p = 1) < 1. (59) 

Imposing the above constraint in its extreme form (i.e. with the equal sign) on (57), we find 
the following relation between / and M 

l ± = y exp(^0 (-127rM 2 ± 2v / 36vr 2 M 4 - 9tt 2 M 2 + 4) , (60) 
3nM + 2 4 V / 

and consequently the following constraint on the maximum value of the disk mass, 

V6 
127 

for I = —0.2204. This is a higher maximum value for the disk mass compared to that 
obtained by Morgan and Morgan [2], 

= ^ « 0.2122. (62) 

In other words, addition of the NUT factor as an extra parameter, not only changes the 
static character of the spacetime (to stationary) but also allows for a higher mass of the 
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Ml™- NUT = — ( 3tt 2 - V9n 4 - U-kA 1 » 0.2427, (61) 

127T V / 
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FIG. 1: I as a function of M by velocity constraint (60). The allowed values of I and M lie inside 
the loop formed by the two curves. 



disk to be acceptable physically. Note that this maximum value of mass is obtained for 
a negative value of I which, by the form of the metric (45), is acceptable. In Fig. 1 the 
diagram of I in terms of of M, given by the constraint (60), is shown and only those pair of 
values (I, M) inside the loop formed by the l + and l_ curves are acceptable physically. For 
a fixed disk mass M the velocity of particles as a function of radius is shown in Fig. 2 for 
different values of I. 

B. Stability and energy conditions 

Having discussed the velocity distribution of the disk particles one could also discuss 
stability of particle orbits under radial perturbations. To do so we use the extended Rayleigh 
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FIG. 2: Velocity of disk particles vs radius for a fixed mass M = 0.15 and different I values 



including I = i.e. Morgan-Morgan disk. For two forbidden values of I the velocity diverges and 



exceeds that of light. 



criterion of stability under radial perturbations given by [23,24], 

f > 0. («, 

where L is the angular momentum of disk particles. For circular geodesies in z = plane 
the specific angular momentum is given by; 

1 d0 

L = 9449 + 9t<t>t = —=(gt<t> + 94>4>-rr)- ( 64 ) 
\j9tt at 

Substituting for ^ from (49) it could be seen that L 2 is an increasing function of p and 
satisfies the above criterion leading to the stability of circular orbits (Fig. 3). Since we 
have already confined the domain of the allowed values of I and M to those respecting the 
velocity of light barrier, it is expected that the same domain of values respect the energy 
conditions as well. This is indeed the case for strong and weak energy conditions as we will 
see below but not for the dominant energy condition. To consider the energy conditions for 
MM-NUT disk space we diagonalize the stress-energy tensor to find the principal pressures 
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FIG. 3: Specific angular momentum as a function of p. 



[23]. The Eigenvalue equation for the stress-energy tensor ; 



T b a e = \C a,b = t,<p. 



constitutes the following eigenvalues 



\ ± = ±(t±Vd), 



where 



(65) 



(66) 



D = (T/ - T^f + 4Jpf and T = T* + T*. (67) 

It could be easily checked that for the discriminant D > the eigenvalues A± are the 
azimuthal pressure P and energy density e respectively [23]. In our case, using the non-zero 
components of the energy-momentum tensor (50), the discriminant, 



D = 14400(1 -p 2 )M 6 e ( " 



fVt.V-aWK [Z 2 e W 2 -" 2 )) - 4M 2 e^( 2 -" 2 ))] ; 



is positive definite and the strong energy condition (SEC), 

SEC = e + P = A+ - A_ = \D\* > 0. 
13 



(68) 



(69) 



is always satisfied. On the other hand the weak energy condition (WEC) and the dominant 
energy condition (DEC) could be written equivalently as the following conditions on the 
magnitude of the gravitomagnetic charge /, 



WEC = e = — A_ > =>■ \l\ < 2M t ^ff> ^) ) (70) 

u 4 + 37rMp 2 



DEC ee e > \P\ =► |/| < 2Mj 2 " 37r ^ eW). (71) 

It could be seen that the most restricted situation happens at the rim of the disk, so, to 
ensure the above conditions all over the disk we apply them at p = 1 i.e., 



WEC = \l\ < 21^/ JT ^exp(— ), (72) 



DEC^l<2MJ 2 -^ exp(^). (73) 
11 ~ V 2 + 3ttM FV 4 ; y J 

In Fig. 4 the DEC, WEC and velocity conditions are shown in a diagram for the same range 
of values of I and M as in Fig. 1. It is noted that the WEC is already satisfied by those 
values of I and M respecting the velocity condition, but DEC could be violated for some 
values of I and M already allowed by the velocity condition and vice versa. We also note 
from (73) that, for the static Morgan-Morgan disk (/ = 0), DEC is equivalent to the velocity 
condition, in other words, the maximum allowed mass is given by M max = as in (62). 

C. Gravitational redshift 

Using the usual formula for the gravitational redshift in an asymptotically flat spacetime, 
for the MM-NUT disk space (which is asymptotically flat) we have, 



, q tt (oo) 1 / ,37rM. 9 I 2 . 371 M , 9 , 

i+z=i/ ty = 77=^( exp< — (p - 2)) +4M^ p <-— <' - 2 ») ■ < 74 » 



AM 2 

It is seen that the largest redshift occurs for an emitter at the centre of the disk ; 



Quipo) 1 ( ,3ttM / 2 3txM 

i + w = \l^ = l exp( ^ } + up exp{ -^ ) ' (75) 



AM 2 

which compared to the Morgan-Morgan disk has the extra /-dependent term in the right 
hand side. In Fig. 5 redshift of photons emitted from the center of the disk in terms of NUT 
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Larger I by v=1 condition 
Smaller I by v=1 condition 
WEC with positive I 
WEC with negative I 
DEC with positive I 
DEC with negative I 



FIG. 4: Dominant and weak energy conditions are shown along with the velocity condition all at 
the rim of the disk. It is notable that there are regions for which the DEC is not satisfied but the 
velocity condition holds and vice versa. 

parameter for two different disk masses are given. We note from the above discussions that 
in general for each value of M there is a range of allowed values of I and indeed in Fig. 5 
the masses and their corresponding NUT parameters are chosen such that all are physically 
acceptable. 

VI. SUMMARY AND DISCUSSION 

In this article using the Ehlers transformation as a tool for finding a stationary spacetime 
from a static spacetime and the 1 + 3 approach to spacetime decomposition which yields 
a physical interpretation of the stationary spacetime in terms of the concepts introduced 
under the general name of gravitomagnetism, we found the stationary spacetime of a finite 
disk dual to the static spacetime of a non-rotating disk (without radial pressure) introduced 
by Morgan and Morgan [2]. The exact form of this metric, which we called MM-NUT disk 
space, is given in oblate ellipsoidal coordinates in (45). The axially symmetric nature of 
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FIG. 5: Redshift of photons emitted from the center of the disk as a function of I for two different 
values of disk mass. 

the spacetime and the form of the gravitomagnetic potential (44), infers a gravitomagnetic 
field with components along the polar directions p and 9. The fact that we find the extra 
parameter of the spacetime in the combination aM suggests that the parameter a could 
be interpreted as the gravitomagnetic monopole charge per unit mass. It is shown that the 
MM-NUT stationary disk admits, through the introduction of the NUT factor, a higher 
maximum disk mass than the static Morgan- Morgan disk for a negative value of I. It is 
found that this spacetime is singular either in the positive or the negative half axis, leading 
to the existence of closed timelike curves (CTCs) in the singular sector. We also showed that 
MM-NUT disk has stable particle orbits and satisfies the strong energy condition. The weak 
energy condition on the other hand holds as long as the values of mass and NUT parameter 
are such that the particle velocities do not exceed that of light. Interesting enough, it is 
shown that the dominant energy condition could be violated for some values of / and M 
which respect the velocity barrier condition and vice versa. Also as a by-product and as an 
example of the application of the Ehlers transformation, we showed that the seed metric for 
both NUT and pure NUT (M = 0) spaces is the Schwarzschild metric. 
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